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ABSTRACT 

The  emerging  theory  of  compressed  sensing  (CS)  provides  a  univer¬ 
sal  signal  detection  approach  for  sparse  signals  at  sub-Nyquist  sam¬ 
pling  rates.  A  small  number  of  random  projection  measurements 
from  the  received  analog  signal  would  suffice  to  provide  salient  in¬ 
formation  for  signal  detection.  However,  the  compressive  measure¬ 
ments  are  not  efficient  at  gathering  signal  energy.  In  this  paper,  a 
set  of  detectors  called  subspace  compressive  detectors  are  proposed 
where  a  more  efficient  detection  scheme  can  be  constructed  by  ex¬ 
ploiting  the  sparsity  model  of  the  underlying  signal.  Furthermore, 
we  show  that  the  signal  sparsity  model  can  be  approximately  es¬ 
timated  using  reconstruction  algorithms  with  very  limited  random 
measurements  on  the  training  signals.  Based  on  the  estimated  signal 
sparsity  model,  an  effective  subspace  random  measurement  matrix 
can  be  designed  for  unknown  signal  detection,  which  significantly 
reduces  the  necessary  number  of  measurements.  The  performance 
of  the  subspace  compressive  detectors  is  analyzed.  Simulation  re¬ 
sults  show  the  effectiveness  of  the  proposed  subspace  compressive 
detectors. 

Index  Terms —  Subspace,  compressed  sensing,  detection 

1.  INTRODUCTION 

Compressed  sensing  provides  a  new  framework  to  jointly  measure 
and  compress  a  sparse  signal  for  sensors  that  need  less  sampling 
resources  than  traditional  approaches.  A  signal  x  £  is  K  sparse 
on  some  basis  'k  =  ^  'ip if  x  can  be  represented  by  a 

linear  combination  of  K  vectors  from  rp  with  K  N.  Given  the 
sparse  signal  x,  the  theory  of  compressed  sensing  shows  that  x  can 
be  recovered  from  M  random  measurements  with  high  probability 
when  M  —  CKXogN  N,  where  C  >  1  is  the  oversampling 
factor  [1],  [2],  The  measurements  are  given  by  y  =  <E>x,  where  <E> 
h  a  M  X  N  i.i.d.  random  projection  matrix  with  each  entry  taken 
from  an  i.i.d.  random  distribution.  Note  that  for  illustrative  purposes, 
we  represent  the  signal  and  the  projection  waveform  in  digital  form. 
However,  it  should  be  clear  that  only  the  projection  results  y  are  to 
be  measured  and  digitalized. 

In  addition  to  signal  reconstmction,  the  CS  framework  can  also 
provide  a  universal  measurement  approach  for  signal  detection  and 
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classification  [3],  [4]  with  a  reduced  set  of  measurements.  A  small 
number  of  random  measurements  suffice  to  provide  sufficient  signal 
statistics  for  the  detection  problem.  The  Matching  Pursuit  algorithm 
is  proposed  in  [3]  to  detect  relevant  signal  components  from  com¬ 
pressive  measurements.  Given  the  compressive  measurements,  the 
compressive  detector  is  proposed  in  [4]  for  signal  detection  and  esti¬ 
mation  without  reconstructing  the  signal. 

The  i.i.d.  random  measurement  scheme  for  compressive  detec¬ 
tion  provides  universality  for  signals  with  different  stmcture.  How¬ 
ever,  since  the  compressive  measurements  are  not  tailored  to  the 
underlying  signals,  the  detector  is  not  efficient  at  gathering  signal 
energy  and  thus  the  performance  is  inferior  to  traditional  detectors 
when  M  <  N.  In  this  paper  we  show  that  the  performance  of 
detection  based  on  compressive  measurements  can  be  significantly 
improved  by  exploiting  the  underlying  signal  stmcture,  leading  to 
the  requirement  of  far  fewer  measurements.  A  subspace  compres¬ 
sive  measurement  matrix  can  be  constmcted  based  on  the  estimated 
signal  subspace  model.  The  analog  projection  waveforms  generated 
through  the  use  of  this  matrix  are  more  efficient  at  gathering  signal 
energy,  in  turn,  leading  to  improved  detection  performance. 

The  subspace  compressive  detector  is  constructed  based  on  the 
theory  of  compressive  sensing,  which  also  draws  on  elements  of  the 
mixed-signal  architecture  introduced  in  [5],  [6],  It  is  different  from 
detection  based  on  classes  of  liner  transforms  [7]  in  that  the  number 
of  measurements  is  flexible  with  more  measurements  leading  to  bet¬ 
ter  detection  performance.  The  subspace  compressive  detector  also 
differs  from  traditional  subspace  detector  [8]  in  that  no  Nyquist  sam¬ 
pling  of  the  received  signal  is  required.  The  analog  input  signal  is 
projected  onto  analog  projection  waveforms  and  only  the  small  set 
of  projected  results  are  sampled  and  digitalized.  It  can  be  shown  that 
far  fewer  samples  are  required  for  the  same  detection  performance. 
Compared  with  compressive  detectors  in  [4],  the  subspace  compres¬ 
sive  detector  achieves  the  same  detection  performance  with  fewer 
measurements  (larger  compression  ratio).  Although  the  subspace 
compressive  detector  does  not  provide  universality  for  signal  detec¬ 
tion  w.r.t  all  signals  in  the  TZ^  space,  it  does  provides  universality 
for  signal  detection  w.r.t  all  signals  in  a  specific  subspace. 

The  sparse  signal  structure  can  be  estimated  with  high  confi¬ 
dence  with  relatively  small  additional  cost.  It  has  been  shown  that 
the  signal  sparsity  model  can  be  estimated  from  limited  i.i.d.  ran¬ 
dom  measurements  of  the  training  signals  [3], [9].  Furthermore,  the 
detection  performance  is  robust  to  the  estimation  error  of  the  signal 
model.  The  Basis  Pursuit  Denoising  algorithm  [10]  is  introduced  in 
this  paper  to  estimate  the  signal  stmcture  for  its  robustness  to  noise. 
The  estimated  signal  model  is  then  used  to  constmct  a  subspace  mea¬ 
surement  matrix  for  subsequent  signal  detection.  The  setting  for  the 
detection  problem  follows  the  typical  CS  setting  where  only  the  ba- 
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sis  'I'  is  available. 


2.3.  Subspace  Compressive  Detector 


2.  DETECTION  OE  KNOWN  SPARSE  SIGNAL 

Here  we  assume  an  explicit  subspace  linear  model  for  the  sparse 
signal  X.  The  K  vectors  of  'P  constmct  a  N  x  K  matrix  H  = 
\ib  ,ib  1,  where  Tii  G  |1,  2, . . . ,  for  i  =  1, . . . ,  if. 

ni  — n2  —'>^K  ‘  j 

The  signal  x  is  represented  as  x  =  HO,  where  0  is  a  if  x  1  vector 
with  all  non-zero  entries.  P/r  =  is  the  projection 

matrix  for  the  subspace  in  which  x  lies.  In  the  following,  we  discuss 
whether  the  detection  performance  can  be  improved  based  upon  a 
priori  knowledge  of  the  signal  sparsity  model  x  =  HO. 

2.1.  Traditional  Detector 

Assume  the  signal  x  G  TZ^  is  known  and  let  n  ~  A/’(0,  cr^Ijv)  be 

1.1. d.  Gaussian  noise.  The  detection  problem  is  to  distinguish  two 
hypothesis  Ho  and  Hi: 


The  measurement  scheme  of  the  compressive  detection  in  Sec.  2.2 
does  not  exploit  the  inherent  stmcture  of  the  sparse  signal.  In  the 
following,  a  subspace  compressive  detector  is  proposed  where  the 
measurement  matrix  <I>  can  be  tailored  to  the  signal  stmcture  so 
that  fewer  measurements  and  better  detection  performance  can  be 
achieved.  Later  we  show  that  the  signal  structure  can  be  estimated 
using  reconstruction  algorithms. 

If  H  and  9  are  known,  the  proposed  subspace  measurement  ma¬ 
trix  is  given  by: 

$  =  G(H^H)"^H'^,  (6) 

where  G  is  an  M  x  K  i.i.d.  random  matrix  with  M  <  K.  The 
detection  problem  is  to  distinguish  between  two  hypothesis  Ho  and 
Hi: 

Ho  :  y  =  $n, 

Hi  :  y  =  $(x-|-n).  (7) 


Ho  ■■  y  =  n, 

:  y  =  x-|-n.  (1) 


It  is  well  known  that  the  optimal  detector  is  the  matched  filter  [11]. 
A  sufficient  statistic  is  given  by:  t  —  (y,  x).  The  performance  of  the 
matched  filter  is  given  by: 


Pd  (a) 


Q 


Q~ 


(2) 


where  Q{x)  =  (27r)“5  e~^  and  Po  is  the  probability  of 

detection.  Here,  a  is  set  such  that  the  probability  of  false  alarm  is 
Pfa  =  a.  The  detection  performance  is  unchanged  even  the  explicit 
signal  model  x  =  HO  is  known. 


It  is  easy  to  show  that  the  sufficient  statistic  is  given  by: 

t  =  =  y^[G(H'^H)"^G'^]“^G0.  (8) 

The  detector  performance  is  given  by: 


Pd  (a)  =  Q 


X^PaX 


(9) 


where 

Furthermore,  we  have: 


x^P§x  = 

=  0'^G^[G(H^H)"^G^]“^G0. 


2.2.  Compressive  Detector 

The  detection  of  transient  signals  that  are  wideband  and  nonstation¬ 
ary  requires  sampling  at  the  Nyquist  rate.  The  sampling  rate  can 
be  reduced  within  the  CS  framework  as  follows.  Given  an  M  x  N 
i.i.d.  random  measurement  matrix  $  with  M  <  N,  where  the  en¬ 
tries  of  <&  are  drawn  from  an  i.i.d.  random  distribution,  the  detection 
problem  with  compressive  measurements  for  a  known  signal  x  is  to 
distinguish  between  two  hypothesis  Ho  and  Hi: 


When  M  =  K,  G  is  invertible  with  probability  1,  which  leads 
to:  x^P§x  =  =  x^x.  Compared  with  (2),  there  is 

no  performance  loss  if  M  —  K.  Note  that  GG^  ^  KImxm, 
G^G  Ri  MIk  xK-  If  'k  is  an  orthogonal  basis,  it  is  easy  to  show 
that:  x^P§x  (M/if)x^x.  Then,  with  M  <  K,  the  detector 
performance  can  be  approximated  as: 


Pd  (a)  ~  Q 


Q-\a)-o/WK\  ^ 


(10) 


Ho  :  y  =  'I’n, 

Hi  :  y  =  ^(x -f  n).  (3) 

The  sufficient  statistic  is  given  by:  i  —  y^($'l>^)“^'l>x.  The  per¬ 
formance  of  the  compressive  detector  is: 


Pd  (a)  =  Q 


'^P<[.X 


(4) 


where 

P<i,  = 

The  detection  performance  can  be  approximated  as  [4]: 


Q 


Q-\a)-^/M/Ni  — 


(5) 


Since  M/N  <  1,  the  compressive  detector  reduces  the  number  of 
measurements  at  the  cost  of  increased  miss  probability.  It  is  also 
clear  that  the  detection  performance  of  the  compressive  detector  can 
not  be  improved  upon  using  a  priori  knowledge  of  the  signal  sparsity 
model. 


The  approximation  of  P\i{a)  by  (10)  is  also  tme  when  ■i/’iS  in  4'  are 
only  approximately  orthogonal  to  each  other.  Compared  with  (5), 
the  proposed  detector  provides  better  detection  performance  with  the 
same  number  of  measurements.  Note  that  only  M  <  K  measure¬ 
ments  are  required.  Much  like  the  compressive  detector  provides  a 
universal  detection  scheme  for  signals  in  the  A^-dimensional  space, 
the  subspace  compressive  detector  exploiting  the  linear  signal  model 
provides  a  universal  detection  scheme  for  all  signals  in  the  same  sub¬ 
space.  The  introduction  of  the  random  matrix  makes  the  detector 
robust  to  magnitude  variations  over  9. 

3.  DETECTION  WITH  UNKNOWN  PARAMETERS 

Very  often  at  the  time  of  detection,  the  signal  is  unknown  or  the  sig¬ 
nal  has  unknown  parameters.  In  the  following,  it  is  assumed  that  the 
signal  sparsity  model  9  is  unknown  but  H  is  available.  The  develop¬ 
ment  of  a  general  likelihood  ratio  test  (GLRT)  detector  with  a  signal 
linear  model  that  samples  the  received  signal  at  the  Nyquist  rate  can 
be  found  in  [11].  In  the  following,  a  detector  that  uses  subspace 
compressive  measurements  is  discussed. 
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3.1.  Compressive  GLRT  Detector 

For  a  compressive  detector  with  i.i.d.  random  measurements,  we 
usually  assume  that  M  >  K,  or  the  detection  performance  is  not 
acceptable.  In  this  case,  the  unknown  parameters  9  in  the  signal 
sparsity  model  can  be  estimated.  With  the  same  hypothesis  testing 
problem  as  in  (3)  and  with  the  a  priori  information  of  H,  it  can  be 
shown  that  the  compressive  GLRT  detector  for  the  unknown  signal 
is  to  decide  Hr  if  the  sufficient  statistic  satisfies: 


when  'I'  is  (approximately)  orthogonal.  Compared  with  (10),  there 
is  performance  loss  due  to  unknown  0  when  x^x/cr^  <  2M.  How¬ 
ever,  with  the  same  number  of  measurements,  the  subspace  com¬ 
pressive  GLRT  detector  outperforms  the  compressive  GLRT  detector 
that  uses  random  compressive  measurements,  as  will  be  illustrated  in 
the  simulations. 

4.  ESTIMATION  OF  SIGNAL  SPARSITY  MODEL 


f'(x)  =  y^C-'(V^C-iV)-'V^y  >  7,  (1 1) 

where  C  =  and  V  =  d?!!.  In  deriving  the  esti¬ 

mate  of  9,  ^  —  (V^C“^V)“^V^C“^y,  is  used.  The  detection 
performance  is  given  by: 

Pfa  =  (7) 

=  '5^)1  (A) 

where  Qx\  the  right-tail  probability  of  a  chi-squared  random  vari¬ 
able  with  K  degrees  of  freedom;  Q  /s.t,  is  the  right-tail  probabil- 

Xk'.A) 

ity  for  a  noncentral  chi-squared  random  variable  with  K  degrees  of 
freedom  and  noncentrality  parameter 

^ - P - (12) 

It  can  be  shown  that  with  a  priori  information  of  the  subspace  (H), 
the  performance  of  the  compressive  GLRT  detector  is  improved 
compared  to  the  case  when  H  is  unknown.  However,  the  improve¬ 
ment  is  limited.  A  more  efficient  detection  scheme  is  needed. 


3.2.  Subspace  Compressive  GLRT  Detector 


As  in  Sec.  2.3,  the  subspace  measurement  matrix  $  is  designed  ac¬ 
cording  to  (6)  and  only  M  <  K  measurements  are  needed.  With  the 
same  hypothesis  testing  problem  (7),  the  likelihood  ratio  is: 


L{y)  = 


p(y;$Hg,7fi) 

piy;'Ho) 


(13) 


where  <I>H0  is  the  MLE  of  under  Hr.  Clearly,  with  M  < 
K,  the  MLE  of  tpHS  is  y.  The  sufficient  statistic  of  the  subspace 
compressive  GLRT  detector  reduces  to: 


t  = 


>  7- 


(14) 


The  detection  performance  is  given  by: 

Pfa  =  Qy2  (7) 


Po  =  Q' 


x^aCA)^’ 


where  the  noncentrality  parameter  A  is  given  by: 


A  = 


G-'  G(H-'  H)-"G' 


-ir^Ti-r 


G9 


(15) 


It  can  be  shown  that  for  x^P,j,x/(t^  ^  M,  the  detection  perfor¬ 
mance  can  be  approximated  as: 


Pn(a) 


Q  ^(«)  - 


Q  ^(«)  - 


[M  x'^x 

T 


(16) 

(17) 


In  the  traditional  CS  setting,  the  signal  x  is  assumed  to  be  K  sparse 
on  a  known  basis  But  the  location  of  the  K  relevant  vectors  of 
a:  in  'k  is  unknown,  which  leads  to  the  matrix  H  being  unavailable 
at  the  time  of  detection.  However,  if  training  signals  can  be  pro¬ 
vided  and  the  compressive  measurements  of  the  training  signals  are 
available,  then  based  on  the  random  measurements  of  these  signals, 
reconstmction  algorithms  such  as  Matching  Pursuit  and  Basis  Pur¬ 
suit  Denosing  can  be  employed  to  estimate  the  locations  of  the  K 
relevant  vectors  in  T'.  Although  these  algorithms  were  initially  pro¬ 
posed  for  sparse  signal  reconstruction,  reconstmction  of  the  sparse 
signal  with  high  precision  is  not  possible  with  very  limited  i.i.d.  ran¬ 
dom  measurements  from  the  training  signals.  Instead,  we  expect  to 
use  these  algorithms  to  get  approximate  information  about  the  signal 
subspace.  Far  fewer  measurements  are  required  for  signal  structure 
estimation  than  for  exact  sparse  signal  reconstmction  [3].  For  the 
purpose  of  detection,  approximate  subspace  information  can  suffi¬ 
ciently  lead  to  great  improvement  on  detection  performance. 

Considering  the  sparse  nature  of  the  underlying  signal  and  the 
wideband  noise  effect,  the  Basis  Pursuit  Denoising  (BPDN)  algo¬ 
rithm  is  introduced  in  this  paper  to  identify  the  sparse  signal  from 
the  noisy  measureemnts  by  solving  the  following  problem: 

mjn  ||'l>'k0  — y„||2-l-A||0||i  subject  to  yn  =  3>('k0-|-n),  (18) 

where  A  <  ||2($'l>)^y||oo  [12],  With  increased  computational  com¬ 
plexity,  BPDN  can  give  better  estimation  for  the  sparse  signal  sub¬ 
space  than  Matching  Pursuit. 

5.  SIMULATIONS 

In  this  section,  the  performance  of  the  proposed  subspace  compres¬ 
sive  detectors  are  evaluated  through  several  simulations.  Only  the 
detection  of  sparse  signals  with  unknown  parameters  are  simulated. 
For  all  the  simulations,  the  sparse  signal  x  is  given  by:  x  =  H0, 
where  H  G  ,  9  G  and  K  <^^N. 

In  the  first  simulation,  the  performance  of  detecting  a  signal  x 
with  unknown  coefficients  9  is  evaluated.  The  SNR  for  the  detection 
is  15  dB.  The  subspace  matrix  H  has  dimensionality  N  —  2048 
and  K  —  200.  It  is  assumed  that  H  is  known  in  the  signal  spar¬ 
sity  model.  Each  entry  of  H  and  9  is  drawn  from  an  i.i.d.  normal 
distribution.  The  column  vectors  in  H  are  approximately  orthogo¬ 
nal  to  each  other.  The  subspace  compressive  GLRT  detector  pre¬ 
sented  in  Sec.  3.2  employs  subspace  random  measurements  that  are 
tailored  to  the  signal  subspace  model  with  the  number  of  measure¬ 
ments  Ml  =  120.  The  compressive  GLRT  detector  presented  in 
Sec.  3.1  uses  i.i.d.  random  measurements  with  the  number  of  mea¬ 
surements  M2  =  400,  and  exploits  the  signal  stmcture  when  es¬ 
timating  the  unknown  coefficients  9.  Simulation  results  in  Fig.  1 
show  that  the  subspace  compressive  GLRT  detector  achieves  better 
detection  performance  with  fewer  measurements.  The  compressive 
detector  exploiting  the  signal  structure  achieves  better  performance 
than  the  compressive  detector  that  does  not  have  the  signal  structure 
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Fig.  1.  Probability  of  detection  of  subspace  compressive  GLRT  de¬ 
tector. 
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Fig.  2.  Performance  of  BPDN  with  subspace  compressive  GLRT 
detector. 
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information  (modeless).  The  simulation  also  shows  that  the  approx¬ 
imations  given  by  (16)  and  (17)  are  quite  accurate. 

In  the  second  simulation,  the  performance  of  the  subspace  com¬ 
pressive  GLRT  detector  that  uses  i.i.d.  random  measurements  to  es¬ 
timate  the  signal  structure  is  investigated.  The  basis  'P  is  a  iV  x  N 
wavelet  Daubechies-4  orthonormal  basis  with  N  —  1024.  The  sig¬ 
nal  is  X  =  H0  where  H  is  composed  of  iT  =  50  randomly  selected 
column  vectors  from  rlr.  Each  entry  of  9  is  drawn  from  an  i.i.d.  nor¬ 
mal  distribution.  The  detection  of  received  signals  with  SNR=  16 
dB  is  simulated.  Signals  are  present  or  absent  with  equal  probabil¬ 
ity.  At  the  time  of  detection,  only  4'  and  K  are  assumed  known. 
For  each  burst  transmission,  200  training  signals  are  first  transmit¬ 
ted.  For  each  training  signal,  the  number  of  measurements  equals 
AK.  The  BPDN  algorithm  is  employed  to  estimate  H  from  averaged 
noisy  measurements.  The  estimated  H  is  then  used  to  constmct  an 
M  X  N  subspace  measurement  matrix  with  M  =  45  according  to 

(6).  With  10, 000  simulation  results  at  each  fixed  Ppa,  the  perfor¬ 
mance  of  the  detector  is  compared  with  the  subspace  compressive 
GLRT  detector  where  H  is  assumed  known.  It  is  shown  in  Fig.  2 
that  BPDN  is  effective  at  estimating  the  signal  stmcture  and  the  per¬ 
formance  loss  due  to  unknown  signal  stmcture  is  small. 

6.  CONCLUSION 

In  this  paper,  we  evaluate  the  perfonnance  of  a  set  of  subspace  com¬ 
pressive  detectors  that  exploit  the  signal  sparsity  model  explicitly. 
FurtheiTnore,  algorithms  are  introduced  to  estimate  the  signal  spar¬ 
sity  model  for  subsequent  unknown  signal  detection.  The  results 
presented  in  this  paper  can  also  be  easily  extended  to  the  case  when 
the  noise  variance  is  unknown. 

7.  REFERENCES 

[1]  E.  Candes,  J.  Romberg,  and  T.  Tao,  “Robust  uncertainty  prin¬ 
ciples:  Exact  signal  reconstruction  from  highly  incomplete  fre¬ 
quency  information,”  IEEE  Trans.  Inf.  Theory,  vol.  52,  no.  2, 
pp.  489-509,  Feb.  2006. 

[2]  D.  Donoho,  “Compressed  sensing,”  IEEE  Trans.  Inf.  Theory, 
vol.  52,  no.  4,  pp.  1289-1306,  Apr.  2006. 


[3]  M.  Duarte,  M.  Davenport,  M.  Wakin,  and  R.  Baraniuk,  “Sparse 
signal  detection  from  incoherent  projections,”  in  Proc.  IEEE 
ICASSP,  2006,  vol.  3,  pp.  305-308. 

[4]  M.  A.  Davenport,  M.  B.  Wakin,  and  R.  G.  Baraniuk,  “Detec¬ 
tion  and  estimation  with  compressive  measurements,”  Tech. 
Rep.,  Dept,  of  ECE,  Rice  University,  2006. 

[5]  S.  Hoyos  and  B.  M.  Sadler,  “Ultra-wideband  analog-to-digital 
conversion  via  signal  expansion,”  IEEE  Trans.  Veh.  Techno!., 
vol.  54,  no.  5,  pp.  1609-1622,  Sep.  2005. 

[6]  S.  Hoyos,  B.  M.  Sadler,  and  G.  R.  Arce,  “Broadband  multi¬ 
carrier  communication  receiver  based  on  analog  to  digital  con¬ 
version  in  the  frequency  domain,”  IEEE  Trans.  Wireless  Com- 
mun.,  vol.  5,  no.  3,  pp.  652-  661,  Mar.  2006. 

[7]  B.  Friedlander  and  B.  Porat,  “Performance  analysis  of  transient 
detectors  based  on  a  class  of  linear  data  transfoims,”  IEEE 
Trans.  Inf.  Theory,  vol.  38,  no.  2,  pp.  665-673,  Mai'.  1992. 

[8]  L.  L.  Scharf  and  B.  Friedlander,  “Matched  subspace  detectors,” 
IEEE  Trans.  Signal  Process.,  vol.  42,  no.  8,  pp.  2146-2157, 
Aug.  1994. 

[9]  Z.  Wang,  G.  R.  Arce, ,  B.  M.  Sadler,  J.  L.  Pai'edes,  and  X.  Ma, 
“Compressed  detection  for  pilot  assisted  ultra-wideband  im¬ 
pulse  radio,”  in  Proc.  IEEE  Int.  Conf  on  Ultra-Wideband,  Sin¬ 
gapore,  Sep.  2007. 

[10]  S.  S.  Chen,  D.  L.  Donoho,  and  M.  A.  Saunders,  “Atomic  de¬ 
composition  by  basis  pursuit,”  SIAM  Set  Comput.,  vol.  20,  no. 
1,  pp.  33-61,  Aug.  1998. 

[11]  S.  M.  Kay,  Fundamentals  of  Statistical  Signal  Processing- 
Detection  Theory,  vol.  2,  Prentice  Hall  PTR,  1998. 

[12]  K.  Koh,  S.  Kim,  and  S.  Boyd,  “11  _ls  :  A  matlab  solver  for  large- 
scale  ll-regulai'ized  least  squares  problems,”  Stanford  Univer¬ 
sity,  Mar.  2007. 


3876 


